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Our Galaxy resides in the center of a vast "Halo" of Dark Matter (DM). This concentration pro- 
duces, in many viable particle physics models, an indirect Weakly Interacting Massive Particle 
(WIMP) annihilation signal that peaks in the Fermi-LAT's energy range. Our knowledge of the 
diffuse background is essential to placing reasonable limits on the DM mass and cross-section. 
We incorporate a systematic variation of the GALPROP galactic diffuse background model, con- 
strained by current cosmic-ray measurements, into a profile likelihood analysis and present pre- 
liminary upper limits on the DM annihilation cross-section using the Fermi-LAT data. 
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1. Signal 

The Milky Way rests at the center of a massive, (^(10^^) [Diemand et al. (2007)], Halo of 
Dark Matter (DM). We have a good sense of the large-scale mass distribution of this Halo from 
N-body simulations, and can approximate that with a simple density profile function. For example, 
here we use an Einasto profile, where p(r) oc exp(— Ar^). If this Halo is comprised of Weakly 
Interacting Massive Particles (WIMPs) capable of annihilation or decay into Standard Model (SM) 
particles, the Eermi-LAT provides a unique all-sky energy window to search for ganmia-ray signal 
resulting from that process. 

To capture the uncertainty in our DM signal, we consider WIMPs of mass 25 to 500 GeV 
annihilating to a variety of SM channels, including bb, ti, t+t~, and . And while N-body 

simulations are powerful, their inability to resolve small scale structure leaves another uncertainty 
in the "boost" that comes from extrapolating their mass functions down to the free-streaming cut- 
off. We therefore use three different mass function exponents (1.8,1.9,2.0) to cover the reasonable 
[Diemand et al. (2008)] range of possibilities. 

While channels like bb deliver their ganmia-ray signal mostly through decay, others (e.g. 
require an extra step in calculating the predicted Eermi-LAT signal. This is the propaga- 
tion and interaction of the final SM decay products. This calculation is performed with the same 
methodology as the rest of our cosmic -ray (CR) backgrounds discussed in the next section. 

2. Diffuse Backgrounds 

Since we expect the Halo signal to be large-scale and slowly varying, masking out regions of 
the sky containing point sources greatly simplifies our analysis without much affecting our sample 
size. We mask the entire galactic plane to ±10°, and Eermi sources with an energy-dependent sized 
mask corresponding to the LAT PSE. 

After masking, we contend only with diffuse backgrounds to the Halo signal: extragalactic, 
instrumental, point source residuals, and galactic diffuse. The first two enter the fit as completely 
isotropic signals with fixed energy spectra [Abdo et al. (2010)]. Point source residuals, products 
of the extended tails of the PSE that reach outside our masking scheme, are modeled using Eermi- 
LAT software ^ tuned to the Eermi Eirst Source Catalogue [Abdo et al. (2010)] measured source 
parameters. 

Originating from interactions (bremsstraulung, inverse compton, and ttq decay) of CR with the 
gas, light, dust, and magnetic fields in the Milky Way, the galactic diffuse is by far the most difficult 
background to model. GALPROP [Strong et al. (2000)] is a code designed to self-consistently 
solve the transportation and interaction of CR within the galaxy from source injection to arrival at 
Earth radius. We assess the quality of a particular realization of GALPROP's parameter space by 
comparing with available CR measurements. This is advantageous in the respect that it makes our 
model independent of the ganmia rays we wish to probe for DM signals. 

The large number and uncertainty of parameters in the GALPROP framework, however, re- 
quires us to make a thorough investigation of that space to quantify our systematics. To that end, 

^ http://Fermi.gsfc.nasa.gov/ssc 
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Parameter 


Range 


Diffusion Coefficient 
Halo Height 
Diffusion Index 
Alfven Velocity 
Electron Injection Index 
Nucleon Injection Index (Low) 
Nucleon Injection Index (High) 
Source Distribution 


1x10^^^4x10^^ 
1 ^ llkpc 
0.33, 0.50 
O^SOkms-^ 
1.8^2.5 
1.7^2.6 
2.26, 2.43 

Parameterized, SNR, Pulsars 



Table 1: GALPROP parameter space. 



we fit with models generated within the parameter space in Table 1. Each model's validity is quan- 
tified by a fit to i^Be/^Be, B/C, and proton data from the HEAO-3, IMP, ATIC-2, CREAM, 
ACE, ISOMAX, AMSOl, CAPRICE, and BESS experiments. The explicit form is, 

;t2 = i,if^i^^^M, (2.1) 

where D and T are model and data respectively, comes from uncertainty in the solar 
modulation (taken to be zhlOOMV), and the sum is over all points and experiments. 

One uncertainty obviously missing from the parameter scan is the normalization of the primary 
electron sources. The sensitivity of a local measurement to nearby sources (or lack of) makes this 
a particularly difficult value to set a parameter range for. Running each model with and without 
primary electrons and finding the difference isolates their contribution so that we can let it float 
freely in our fits to the data. 

3. Profile Likelihood 

After binning into 12 annular and 80 logarithmic energy bins, we compare our model with the 
Fermi-LAT gamma-ray data and find the Poisson Likelihood, L [Mattox et al. (1996)]. For each 
DM channel and mass we then produce a Maximum Likelihood, 

LjiOuM) = P^^j(^/;S:^flX5 0DM), (3.1) 

i 

by selecting the best- fitting linear parameters, a (EGB, instrumental, and primary electron 
normalizations), for each DM normalization, Odm- The j represents the diffuse model and the 
product is over each bin, /. 

Finally, we convert L into a Test Statistic for every GALPROP model that includes the CR-;^^ 
(Eq. 2.1) information, 

Pf%{0DM) 
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Figure 1: Qualitative plot of an upper limit set using the Profile Likelihood method on the Test Statistic 
defined in Eq. 3.2. Each curve represents a different GALPROP model. 




Figure 2: A qualitative plot of an under-sampled minimum in the Profile Likelihood. 



This leaves us with a family of curves (Fig. 1), the minimum contour of which comprises the 
Profile Likelihood. The contour should behave as a with one degree of freedom, and we set our 
95% confidence upper limit on 0dm to where it rises above the minimum by 3.84. 

To be confident with the limit it is necessary to have sufficiently populated the model space 
such that the minimum is the true minimum, and that the nearby sampling is dense enough to 
smoothly detect the profile out to the upper limit. Under- sampling the parameter space leads to a 
situation like in Figure 2, where the limit is derived from a single model, and both the location of 
the true minimum and the shape of the curve around it are unclear. 

After our initial, completely random, sampling of the parameter space in Table 1, we found 
ourselves in precisely that situation. Figure 3 illustrates the initial sampling as a scatter plot for the 
CR and gamma rays L. The minimum will be near the lower left hand corner, and sampling should 
be densest around it. Clearly this was not the case. 

3.1 Classification Tree Sampfing 

Being limited mainly by computation time, the best way to populate a specific region is to have 
a rough idea what its L will be beforehand. Predicting the gamma-ray fits would be problematic, 
due to the variety of DM signals being explored. The CR-;^^ is much more straight forward, and 
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Figure 3: A scatter plot of the CR vs. gamma ray maximum Likelihood for our present models. 
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Figure 4: Present 95% upper limits on pure bb annihilation from boosted Einasto DM. Also shown are 3(7 
limits from an independent Fermi Collaboration analysis with an unboosted NEW profile, 4 kpc Halo height, 
and SNR-based CR source distribution [Zaharijas et al. (2010)] [Cuoco et al. (2010)]. 



so we use it to train a classification tree. The tree can then be used to pre-select GALPROP models 
that agree with CR data. This will reduce much of the y-axis scatter in Fig. 3, hopefully allowing 
us to sufficiently sample the global minimum with our available computer resources. 

4. Preliminary Results 

Figure 4 shows our present 95% confidence upper limits on a pure bb channel WIIMP. As 
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mentioned earlier, these limits are based on a profile likelihood minimum composed of a single 
GALPROP model. This makes the limits very preliminary and it is unclear in which direction the 
final ones will move. Also shown in Fig. 4 are the closely -related limits derived by an independent 
Fermi Collaboration analysis which uses a single "conservative" GALPROP model with all of its 
gas (target) map contributions free to float. 

A further caveat, there are persistent residuals even in our best fits. We hope to account for 
much of this by correcting known deficiencies to our model, foremost being the isotropic approx- 
imation to the inverse compton signal. We also plan on addressing uncertainty in the gas distribu- 
tions by allowing some float in the collective contributions produced by interaction with molecular, 
atomic, and ionized Hydrogen. Given the level of uncertainties discussed above, these upper limits 
should be considered as purely illustrative of the Fermi LAT potential to constrain DM annihilation 
or decay from the Galactic halo. 

5. Summary and Outlook 

We have conducted a preliminary investigation of possible DM annihilation signal in the Milky 
Way Host Halo, and within the systematic uncertainties of our background model, find no sig- 
nificant detection. We model all known backgrounds, and use the Profile Likelihood method of 
marginalizing over the many nuisance parameters that come from uncertainty in the CR-induced 
galactic diffuse emission. Sampling the model parameter space adequately is a computationally 
difficult task, one we are addressing by intelligently pre-selecting our models with a classification 
tree. We expect the combination of increased sampling density of our maximum likelihood region 
and upcoming improvements to our models (anisotropic inverse compton & inclusion of target un- 
certainty) to yield a robust limit that takes into account the large systematics inherent to a physical 
model of the galactic background. 
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